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SUMMARY 



The investigation described in this report was made 
to study the velocity distribution in an open, in a par- 
tially open, and in a partially expanding jet. The open- 
jet observations reveal minor systematic discrepancies 
from Tollmien's theoretical velocity distribution. The 
shearing- st r ess distribution for the partially open jet 
was determined. The value derived for the ratio of mixing 
distance to jet width was found to be in close agreement 
with the corresponding value for the open-jet boundary. 
The streamline pattern in a partially expanding channel 
was obtained from the observed velocity distribution and 
plotted, after which the distribution of the mixing dis- 
tance for this case was also ascertained. 



INTRODUCTION 



Based on professor Prandtl T s equation for the appar- 
ent shearing stress induced by turbulent momentum inter- 
change (reference l) 
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(where u = mean velocity in x direction; y = velocity at 
right angles to it; p = density; and I = mixing distance), 
W. Tollmien (reference 2) has analyzed various cases of 
free turbulence - i.e., of flows without boundary walls. 
How the purpose of the present investigation was to experi- 
mentally check one of his examples: the uniplanar jet ex- 
pansion upon emergence from a linear slot, and to prove 
the proportionality constant c, which he introduced as 



I = C X 



* u Uber turbul ent e St rahlausbr eitung . M 
vol. V, no. 1, ifS4 t pp. ,42-54. 



(2) 
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in his article. (The corresponding laminar problem has 
"been theoretically analyzed by H. Schlichting, in Zeit- 
schrift fur angewandte Mathematik und Mechanik, No. 13, 

jfr 1933, p. 260; see also W. Muller, in the same issue, p. 

7 395.) In addition, the case of turbulent expansion in a 
^« partially open and in a partially expanding jet, were to 
he checked by the same methods* 

EXPERIMENTAL ARRANGEMENT 



Tollmien's example treats the uniplanar jet expansion, 
to wit: an air stream emerging through an assumedly linear 
slot intermingling with still air. The assumption of lin- 
ear slot is, of course, not exactly realizable in an exper- 
iment, although a slot having a height equal to twenty 
times its breadth, approaches this condition fairly closely 
(fig. l) m A simple calculation then affords a linear slot 
in extrapolation. To render the problem amenable to uni- 
planar treatment, it was necessary to provide such flow 
conditions, for which the velocity distribution was the 
same in all horizontal sections (as far as the upper and 
lower limit range). Compliance with this premise was de- 
manded from the second and third cases also. For the 
first case - that is, perfectly free jet expansion - the 
two-dimensionality was summarily predetermined, whereas in 
the other two cases, it was more difficult. However, be- 
fore going into details, it is advisable to describe the 
experimental set-up itself. 

The air is supplied by a two-stage blower. A honey- 
comb is mounted at G, while aft of it the circular sec- 
tion merges into a rectangular pressure chamber which 
gradually tapers to a slot S. This part of the equipment 
served in all experiments, although in the second and 
third cases some other parts were added. A lateral orifice 
on the pressure chamber D connects with a manometer, 
which is held constant by continuous regulation of the 
r.p.nu of the motor, thus insuring a constant air volume 
in the pressure chamber as well as through the slot. The 
actual test section - after passage of slot - was, in the 
first case, the free air space before the slot. Thus, the 
first problem consisted in dividing this space into exact 
coordinates, so as to be able to measure it point by point. 
This was accomplished with the support shown in figure 2. 
A small table mounted on two parallel posts can be raised 
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or lowered* A horizontal balance bar attached to the lit- 
tle table, supports, on its end, the pressure survey appa- 
ratus intended- to indicate the pressure at the__ part iculafc 
point of the air space. The survey apparatus connects, by 
means of .a pipe running through the balance bar, with a 
rubber tubing which, fastened at the other end of the bar, 
leads to the manometer. The two uprights are mounted on a 
heavy metal support which rests on four setscrews on a 
wooden base. The latter rests only with the four corners 
on the floor, whereby one corner itself is replaced by a 
setscrew to allow for the unevenness of the floor. The 
setscrews of the metal base assure perfect perpendiculari- 
ty of the uprights. The bar carrying the pressure survey 
•device, a pitot tube, was at first set into severe oscilla- 
tions by the air stream, but these could be almost com- 
pletely removed, after various attempts, by the damping L. 
This is effected throiigh a loose piece of metal placed on 
the bar linked to an inclined rocker at a fixed point and 
which slides back and forth with friction when the bar os- 
cillates. The origin of the coordinate system of the test 
range was placed in the intersection of the diagonals of 
the slot for reasons of symmetry of the whole test field - 
direction x in flow direction, direction y, horizontal and 
perpendicular to it. In direction z (vertical), the pres- 
sure grading was, as stated before, constant. 

The total head was recorded by pitot tube and Prandtl 
type manometer. The static pressure was recorded with a 
tube having lateral orifices and an inclined manometer, 
Prom the resultant dynamic pressure, the velocity distri- 
bution was computed after correcting for air pressure and 
temperature. In the second test the creation of a two- 
dimensional flow was, of course, extremely difficult. The 
original shape of the channel was as shown in figure 3. 
But, .-this arrangement proved unsuitable, because the 
suction of the still air, conditioned by the air stream, 
forced it to be sucked around the upper and lower edges 
and the thus-created vortices at the edges soon destroyed 
the uniformity of flow in the vertical. 

This objectionable feature was removed by fitting the 
curved pieces of metal shown in figure 4. Even so , there 
remained difficulties despite the careful alinement of the 
boundary walls (by means of a light beam), until the slot 
was narrowed from 6.5 to 3 cm (2.56 to 1,18 in.) for a 
height of 65 cm (25.6 in.), 'and a second close-meshed screen 
was inserted between pressure chamber and slot. This final- 
ly insured a satisfactory velocity distribution in the ver- 
tical (fig. 5). 
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For the first and second experiments., it was very de- 
sirable to "be able to adjust "both pitot and static tube in 
the pertinent flow direction at the extreme mixing bound- 
ary between jet and inducted air, where the flow direction 
deflects considerably from that of the principal flow f.x 
direction)* The tube had to be made to pivot about the 
test point. The necessary device as shown in figure 6, 
includes a balance bar with a scale division at its outer 
end* Then the static tube, which is very responsive to 
directional changes, was rotated until it indicated maxi- 
mum pressure, after which the corresponding angular setting 
was recorded and the pitot tube pointed in the same direc- 
tion for the total-head measurement. After thus solving 
most of the experimental problems, there still remained 
the difficult feat of attaining a satisfactory static pres- 
sure reading. 

Static Pressure 

The static pressure was recorded, as previously stat- 
ed, with a pitot tube having four symmetrical orifices on 
the sides (fig. ?) . It was soon found that this tube re- 
corded a very marked negative pressure, amounting up to 15 
percent of the dynamic pressure in the main mixing zone of 
air stream and. still air., as against one-half percent and, 
at that, indicative of positive pressure, according to 
Tollmien's theory. Consequently, these measurements, rath- 
er than representing real values, imply that survey devices 
of this type do not record the real static pressure in tur- 
bulent flow. They are extremely responsive to directional 
changes, and even a change of a few degrees causes a pres- 
sure drop in the manometer, which is no longer in accord 
with actual pressure conditions in the jet. 

Since, in turbulent flow, there is a continuous change 
in flow direction, of the individual fluid particles to- 
ward the main flow direction, the manometer reading is 
lower than the true static pressure. And this figure is 
so much less as the turbulent interchange is greater - 
hor.ee the premise, confirmed by personal observational 
study, appears justified. That is, a survey apparatus of 
this kind indicates a minimum of the pressure at the points 
of maximum interchange and consequently, gives the point of 
maximum shearing stress, but for the rest, is totally un- 
suitable for recording the actual static pressure in tur- 
bulent flow (reference 3). 

The next step was to .design a device which was not re- 
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sponsive to changes of flow direction. After various 
abortive attempts, the form shown in figure 8 was adopted. 
The static pressure tap was placed within a • turned-up hol- 
low cylinder, whose walls were hollow and served as pres- 
sure leads to the manometer. A calibration for direction- 
al change revealed it to "be unaffected by flow direction 
over a range of approximately 80°; that is, the reading 
was the same after being turned through a ±40° angle rela- 
tive to the main flow direction (fig. 9). 

On the other hand, it was still markedly influenced 
by the dynamic pressure (fig. 10) which, of course, de- 
creases as the air speed increases. Since most measure- 
ments were made at such speeds, when the effect of dynamic 
pressure is small, the order of magnitude of the static 
pressure was obtainable to such a degree as to make the 
static pressure effect amenable to estimation when comput- 
ing the shearing stresses. 

The interaction of the static pressure on the shear- 
ing stress calculation could be disregarded in first ap- 
proach (as proved elsewhere), hence it afforded a more ex- 
act determination of the static pressure in our particular 
case. In any case, it appears to be desirable, however, 
to concentrate on the design of a survey apparatus which, 
totally unaffected by main flow direction and amount of 
dynamic pressure, affords an unobjectionable static pres- 
sure reading* 

RESULTS OE TESTS 
Op en J et 

The thus-obtained velocity distributions for various 
slot spacings a are shown in figure 11. These profiles 
are plotted nondimens ionally in figure 12 after dividing 
the pertinent velocities by the velocity maximum of the 
particular profile and the related position coordinate by 
the coordinate of the halved maximum velocity. 

The velocity profiles may be expressed by 

Comparing two profiles each for the constants p and 
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q, we have: 

p = - 0.50, q ' = 0.997 

in accord with the theory which stipulates p » - \% q = 
+ 1. 

For comparison with Tollmien's theoretical velocity 
dist ri "but ion , we- included his theoretical curve in figure 
12. . Soth the theoretical and experimental curves have 
"been normalized so as to "be agreeable in the points y = 0, 

j „. / maxN „ _ n max 

* - %ax and y = y v 2 >' • 

Approaching maximum velocity, the measured velocity 

is seen to be higher, and at the jet boundary lower, than 

the theoretical. The same phenomenon was observed in the 

rotationally symmetrical wake (unpublished measurements by 

Hikuradse and Schlicht ing) • The proportionality constant 

c in the mixing distance equation I = c x was defined* 

at c = 0.0165, which gives 1 = 0.0S5, when 2b - total 

b 

jet width* 

Partially Open Jet 

The obtained velocity profiles for different slot dis- 
tances a are shown in figure 13. Whereas the first three 
profiles (a = 10, 20, and 35 cm (3«94, 7.87, and 13.78 in.)) 
still show a slot-width effect and therefore do not count 
in the final evaluation, the others, farther away, do not 
reveal it. At the same time it was important to know 
whether these profiles were mutually similar and the flow 
thus followed a simple law within the entire range. The 
procedure was the same as for the open jet. The individu- 
ally recorded velocity was divided by its momentary maxi- 
mum of the profile, and the individual position coordinates 
by the coordinate of the velocity which, on the free side 
of the profile, equals half the velocity maximum. The 
thus-obtained nondimens ional values for the individual 
profiles were plotted on one sheet. It was found that all 
profiles are coincident in one curve (fig. 14), thus prov- 
ing that the flow in this case also follows a simple law 
of similitude within the entire range. 



♦Owing to the reduction of the finite slot width to linear 
slot, we have x = a + e (a = distance from slot, e = - 
4.5 cm (17.7 in.) • 
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The comparison of two adjacent profiles gave: 

consequently, p = - ^ , q = 0.99 

Hereby it v/as necessary to so reduce the distance 
from the slot of finite width to a linear slot, that a 
constant value e was added to distance a, making the 
theoretical distance x = a + e» The values for p, q t 
and e were defined as follows: ' Since the individual 
profiles are coincident, the momentary velocity must he 
dependent on the theoretical slot' distance as well as on 
a function which- refers the jet width to this distance. 
That is, 




holds again, with 

x = a + e 

a "being the measured distance from the slot, e the 
above-cited constant additive length. Now, since the max- 
imum velocity is simply a function of the slot distance, 
the equation 

(a x + e)^ _ u m (a x ) 
(a a + e)P u m (a a ) 

must hold; whereby a x and a 2 are the measured distance 
of profiles 1 and 2 from the slot, and u m (a x ) and 
u m (a 2 ) , the maximum velocities of profiles 1 and 2. 
Likewise, it must he 

( a 2 + e) p _ 1% (a 8 ) 
(a 3 + e) p u m (a 3 ) 

for profiles 2 and 3. 

7 Prom these equations, e and 
at 

e = 20.1 cm (7.9 in.) , 



p are then computed 

p = - * 
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q may also be expressed by an equation, for according 
to the law of similitude for velocity distribution, equal 

?- values have equal — values. The selection of 



u. 



points in two profiles, so that — = — ■ stipulates 



y (a x ) (a 1 + e) 
y (a 2 ) (a 2 + e) 



when y (a x ) and y (a 2 ) '. denote the y coordinates of 

points ~- h. From this q may be computed. Even so, 
"in ^ 

the obtained value cannot be considered final because the 
shape of the velocity profile is dependent on the wall 
friction as well as on the mixing with the still air. Both 
processes must be treated separately. Substituting for 
this y, measured with inclusion of the friction zone, the 
distance of that point on the free side where the velocity 

equals | -u max and j "^max » sa y> at points 1 and 2 , re- 
moves the wall-friction effect. In this manner the value 
q x = 0.997 was obtained. 

A logarithmic plot of the velocities versus y for 
wall proximity, reveals a rise of the velocity with the 
l/7 power of the wall distance, while for this zone q 2 = 

0.8995. 

This value for q 2 can also be deduced theoretically. 
It has been shown that 



u = c y 



1/7 (3) 



is valid. On the other hand, v. Karman (reference 4) de- 
fined the shearing stress at the wall with 

1/4 7/4 

T o ** P v c = proportional, (4) 

V = kinematic viscosity) 

The momentum theory (equation (9)) stipulates 



0 
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Then it is 



| u = x p .f(n) , where r\ = ~ (6) 

x^- 



that is, with y = x 9 - T) 



x 2 p f 2 (n) x* 



i or 



J Ik ~ x 2p+ q -l ( i 

? P 

To render (3) agreeable with (6) it must be 
f o'r wall proximity, 



or c ~ x' 



which, written in (4) gives: 



1 7 

- ~ u 

P 

From (7) and (8) then follows 



-O-^x 4 ^" 7 ' (8) 



x 2 P+^-l „ x 4 ST 7/' 
or 

/_ 7 \ / IN 

1 = 0 



that is, 



H 5 5 



The insertion of the ascertained value p - - J, then ac- 
tually affords 



q = .q 2 = 0.9 
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Calculation of Shearing Stress- • 
According to the first equation of motion, it is 

3u Su 1 BP 1 3t 

U ~ + V ~ = - — 5— + — 3— 

3x oy ■ p ox p ay 

Integration, according to y gives: 

y . £j _ ^ - y 



T 

p 



which, because g~ = - "becomes 



T 

P 



= 2 / u || dy + vu 



o o 
which may also he written: 



1 r 



7 + L / §E dy 



T 
P 



7 = a / u 2 dy - * ^ / udy + i / g dy 
0 dx o dx 0 p o dx 



(9) 



The differentiation according to x may, on the hasis of 
the law of similitude for velocity distribution 

u = xP tit)) 

y 

"be effected with r\ = -~ • 

X* 

Disregarding the pressure term which, according to 
the measurements, is so small as to he negligible for the 
shearing stress calculation, affords 



T 

P 



7 „ JL [x 2p / f a d 3 _ X V f 3L xP / fdy 
3x 0 ox o 



= (2p+q) x 2 P+<l-l/ f 2 dn - q Tl x 2 P + ^ 1 

- xP(p+q) xP+^l" 1 f /fdn + q x 2 ?**- 1 f S |o- 

0 

= (2p+q) x2p-l /f a dy - (p+q) x 2 P-! f / fdy 
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Thus the shearing stress formula for any value of 
q "becomes: 



T 
P 



C(2p + q) I u 2 dy - (p + q) u / udy] 



and 



do) 



With y m as coordinate of the velocity maximum Vi m , 
gives for the wall region 0 ^ y = y m with p = - -g- and 
q = 0.9 

y y 

(11) 



T 

P 



7 = ~ lf^ 2 d.y + 4u / udy] 

o^o o- 

and for the outside zone y = y with p = - -Ar, q = 1.0 



T 

P 



== - gr- u / udy 
<j x 
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( Ha) 



m 



From (11) and ( 11a) follows the shearing stress in the out- 
side zone at 



T 

P 



y n 1 y m y m 

= - [ / u 2 dy + 4u m / udy + 5u / udy] (12) 

o x o o y m 

Equations (ll) and (12), containing only known values, 
served as "basis of the shearing stress calculation. 

The values of the integrals at the part icular points 
y were graphically integrated. The resulting T/p val- 
ues give the difference in shearing stress at any point t>f 
the flow from the corresponding wall figure which is not 
zero; although the shearing stress at the jet edge rela- 
tive to the free space must "become zero* Plotting the 
values accordingly reveals the shearing stress as passing 
through "kero at the point of maximum flow velocity, as 
equally follows from equation (l2) # This is readily under- 
stood because at this point the shear due to wall friction, 
as well as that due to mixing, must disappear • In figure 
15 the shear distribution is shown plotted against jet 
width. 

Calculation of Mixing Distance^ 

Following the determination of the shearing stress, 
the mixing distance was computed according -to 



p ■ v 



9y 



3u 
3y 
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2"S. values were obtained "by graphical differentiation 

of the curves established for the velocity distribution in 
y direction. The value for the partially open jet is 

1 = 0.0684 

D * 

Here "b = width of mixing zone - that is, exclusive of fric- 
tion zone on the fixed wall. This figure is fairly agreea- 
ble with that for the free jet "boundary = 0.0682^, 
given in Tollmien^ cited report, 

Plow in Partially Expanding Jet 

Experimentally, this case is not unlike that of the 
others except for the closing of the fourth side of the 
previous three-sided channel. The balance "bar carrying^ 
the pressure survey apparatus was inserted through openings 
into the newly installed side wall of the channel. The 
openings for different slot spacings harmonized exactly 
with the profile of the "balance bar. The unused openings 
were carefully plugged. The measurements included the to- 
tal head and the static pressure in eight sections, afford- 
ing the velocity distribution in the channel for the dif- 
ferent slot spacings* (fig. 16). Individual graphical in- 
tegration of the velocity profiles then gave the volume of 
flow for each. The points of equal flow volume were then 
connected, resulting in the flow pattern shown in figure 
17. 

Lastly, we computed the mixing distance for the case 
of abruptly expanding ch a nnel. The velocity profiles be- 
ing no longer similar, the calculation of the mixing dis- 
tance had to be made on the basis of the general equation 
(9). The shearing stress T Q at the wall was computed 



*The measured velocity profiles reveal "bulges 11 near the 
points where the prefix of the velocity changes. This is 
due to the fact that at the point where the turbulent ve- 
locity fluctuation is of the same order of magnitude as 
the mean velocity, the determination of the velocity from 
pitot-tube readings is considerably falsified; that is, 
the velocity is always too high. For the subsequent eval- 
uations, this error was minimized by plotting the velocity 
profiles smooth in the vicinity of the zero points. (Com- 
pare the dashed lines in fig. 16.) 
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from the \jl power law for velocity distribution at 



T, = 0.022 5 P -f- 



1/4 



The integral in (9) was defined graphically from the 
measured velocity distributions. Thus . figure 18 illus- 
trates the obtained distribution of the mixing distance 
over the section for five velocity profiles at distances 
a = 10 to a = 60 cm (a = 3.94 to 23.62 in.) from the 
discharge opening* The notable fact here is that in the 
vicinity of the wall the mixing distance increases consid- 
erably more with the wall distance than in the usual chan- 
nel or pipe flow, where I = 0.4 y (reference 6) in wall 
proximity. It does therefore seem of importance as far as 
the mechanism of turbulence is concerned, whether the air 
sucked in and carried along by the Jet was at rest as in 
both the open and in the partially open jet, or whether 
the air was already in turbulent motion beforehand, as was 
the case here. At greater wall distance the mixing dis- 
tance assumes a constant and practically equal value for 
all profiles; that is, I Z 1.3 cm (0V51 in.) and, since 
b = 18 cm (7.09 in.) , 

A = 0.072 



b = total 
pipe is 

the same. 



channel breadth. '(The 
'IlS = 0.07 for pipe 



value for the 
center - that 



straight 
is, almost 



This method of computing the mixing distance is no 
longer permissible for profiles at greater distances from 
the exit opening because the simple mixing distance for- 
mula (l) no longer holds, as a result of the considerable 
directional differences of the individual streamlines. As 
for the rest, the determination of the mixing distance by 
this method can lay no great claim to accuracy, because 
the flow was not strictly two-dimensional, and the flow 
volume therefore was different in the individual sections ■ 
a fact which is felt particularly disturbing in the graph- 
ical differentiation, according to x. Even so, it may be 
seen that the mixing distance in the suddenly expanding 
channel is of the same order of magnitude as for the open 
and the partially open jet # 
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In conclusion, the writer wishes to acknowledge his 
indebtedness to Professor Prandtl for suggesting the in- 
vestigation, and for- his valuable assistance; and to Dr. 
W. Tollmien and Dr. H. Schlichting, for their critical 
reading of the manuscript. 

Translation by J, Vanier,-* 
National Advisory Committee 
for Aeronautics. 
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Figs. l t 2 f 4 




Figure 1.- Experimental layout 

for open jet measure- 
ment s. (Blower with slot-shaped 
contraction. ) 




Figure 4.- Experimental setup for 
partially open jet. 
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Figure 3.- Partially open 
jet. 




Figs. 3,5,6,7,8,9,10,11 
~t 4— = 



D 



Figure 7.- Static tube 
with side 

orifices. 
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Figure 5.- Velocity distribu- 
tion at right angles 
to main flow plane for partial- 
ly open Jet. 



Figure 8.- Static pressure tube 
with inside holes. 




Figure 6.- Device for rotating the 
pressure survey appara- 
tus about the vertical axis. 
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rigB. 12,13,14,15 
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Figure 12.- Non-dimensional velocity 

distribution in two-di- 
mensional open jet. 
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Figure 14.- Non-dimensional velocity 



Figure 13.- Velocity distribution I 

in partially open jet » 

for different distances from the y 

slot. 1 



Figure 15.- Shear distribution in 
partially open jet. 
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Figs. 16,17,18 





Figure 17.- Flow picture of partially expanding jet. 




Figure 18.- Distribution of mixing distance in 
partially expanding channel . 



